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Abstract 

QhI We consider a bipartite mean-field model in which both the interaction constant and 



the external field take different values only depending on the groups particles belong 
to. We compute the exact value of the thermodynamic limit of the model exploiting 
a tail estimation on the number of configurations that share the same value of 
^ , the magnetization and we analyze the critical points of the pressure functional 

o 

^SJ ' associated to the symmetric version of the model as the external field is away or 

^ : 

Tlj- , small. 



o 

(N 
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Introduction 

Bipartite mean-field model have been introduced since the 50s to reproduce the phase 
transition of some materials called metamagnets. In particular in [6] , [2] and [7j a bipartite 
mean-field model is used to approximate a two-sublattice with nearest neighbor and next- 
nearest neighbor exchanging interactions. The same model has also been used to study 
the loss of gibbsianess for a system that evolves according to a Glauber dynamics [8]. 

In recent times the general version of these models have been proposed in the attempt 
to describe the large scale behaviour of some socio-economic systems [3], assuming that 
individuals decisions depends upon the decisions of others. 
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The investigation of the model introduced in [3J has been pursued at a mathematical 
level in [5]. It has been shown the existence of the thermodynamic limit of the pressure 
exploiting a monotonicity condition on the Gibbs state of the Hamiltonian (see pQ). The 
factorization of the correlation functions has been proved for almost every choice of pa- 
rameters and the exact solution of the thermodynamic limit is computed whenever the 
Hamiltonian is a convex function of the magnetizations. 

In this paper we carry on the analysis of the mathematical property of the bipartite 
mean field model. Firstly we compute the exact solution of the thermodynamic limit. We 
exploit a tail estimation on the number of configurations that share the same value of the 
vector of the magnetization. This technique is the same used by Talagrand to compute 
the thermodynamic limit for the Curie- Weiss model [9J. Then we analyze the critical 
points of the pressure functional associated to a symmetric bipartite mean-field model in 
the case in which the external field is absent or small. This analysis highlights for which 
values of parameters the model undergoes a phase transition. 

This work is organized as follows. Chapter one introduces the notations and states 
the main results. Chapter 2 contains the proofs. The appendix contains the proof of the 
lemmas that make the paper self contained. 

1 Definition and Statement 

We consider a system of particles that can be divided into 2 subsets Pi,P2 with 
Pi n P2 = 0, and sizes \Pi\ = Ni, where Ni + N2 = N. Particles interact with each other 
and with an external field according to the mean field Hamiltonian: 



The (Tj represents the spin of the particle i, while Jij is the parameter that tunes the 
mutual interaction between the particle i and the particle j and takes values according 
to the following symmetric matrix: 




(1) 



2 



Ni N2 



No I 



Jii 



T* 
•J 1 



12 



I v 



112 



>22 



where each block 3 is has constant elements Jis. For / = s, J;; is a square matrix, whereas 
the matrix J^^ is rectangular and all entries can be either positive or negative allowing both 
ferromagnetic and anti-ferromagnetic interactions. The vector field takes also different 
values depending on the subset the particles belong to as specified by: 



N2 < 

. V J 

where each h; is a vector of constant elements hi. The joint distribution of a spin config- 
uration cr = ((Ti, . . . , (Tat) is given by the Boltzmann-Gibbs measure: 

exp(-/3i^Af(cr)) 



PN,J,h{o-} 

where Z7v(J,h) is the partition function 



(2) 



Z^(J,h) = ^exp(-/3if^(o-)). 

cr 

By introducing the magnetization of a set of spins A as: 



mA[cr) 



and indicating by mi{a) the magnetization of the set P;, and by ai = NjN the relative 
size of the set P;, we may easily express the Hamiltonian ([T]) as: 



(3) 



where the function g: 

g{mi{a),m2{cr)) = ^ ( Jiimi(cr)^ + 2aia2Ji2'mi{a)m2{cr) + J22"^2(cr)^ ) + 



2 

+ ai/iimi(cr) + a2^2'^2(c) (4) 
depends from the reduced interaction matrix J and external field vector h: 

J 



Ju 


J121 








W12 


^22^1 







The existence of the thermodynamic limit of the pressure: 

p^(J,h) = llnZjv(J,h) (5) 

associated to the model defined by Hamiltonian and distribution ([2]) is shown in the 
paper [S] exploiting an existence theorem provided for mean field model in In this 
paper we deal with the exact calculation of this limit for all the values of the parameters. 

Proposition 1. Consider the Hamiltonian defined in (J\) and the function: 

/(/ii, /is) = (3g{iJi, H2) - «! J^(/ii) - a2-'^(/U2) (6) 

where g is given by ^ and: 

j^(x) = ^ ((1 + x) ln(l + x) + (1 - x) ln(l - x)) . (7) 

Given parameters Ju, J12, J22, hi, h2 and a, the limit for large N of the pressure PAr(J, h) 
associated to the Hamiltonian and defined in ^ is the following: 

lim pjv(J,h) =ln2+ max f[^i,^2)- 

Therefore to compute the exact value of the thermodynamic limit of the pressure 
PAr(J,h), we have to maximize the function / with respect to the variables /ii and fi2- 
Differentiating / we obtain: 

-(yUi,/i2) = l3{a\Jii^ii + ala2Jl2^^2 + aihi) - ai tanh"-^(/ii) 



(/ii,/i2) = /3(aia2-^i2/ii + celJ22l^2 + Ci2h2) — a;2tanh" (/i2) 



thus the mean field equations of the model are: 

{/ii = tanh ^/3(aiJii /ii + 0:2^12 yU2 + hi) 
/i2 = tanh [(3{aiJi2 yUi + 0:2^22 /i2 + ^2) 
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We observe that considering the symmetric case of the model, in which a = 1/2 and 
Ju = J22, and changing the sign of both the parameter of mutual interaction J12 and 
the first component of the external field hi, we obtain a system of mean field equations 
symmetric to the previous one with respect to the vertical axis. Therefore to analyze the 
critical points of the pressure functional / in the symmetric case, it suffices to consider 
only one sign for the parameter J12. In particular we will study the cases in which Ju and 
J12 have opposite signs in order to keep the graphical approach introduced in [2], useful 
to understand the analytical results. 

Proposition 2. Consider the Hamiltonian defined in (QP whit a = 1/2, Ju = J22 and 
hi = h2 = 0. Denoted with Xm and Xm respectively the bigger and the smaller eigenvalue 
of the reduced interaction matrix J associated to the Hamiltonian, as Ju and J12 have 
opposite signs, the pressure functional f given by ^ admits the following structure of 
critical points: 

1. if Ju > and J12 < (Xm = Ju - Ju] Am = Ju + J12) 

(a) < (3 < y- : the unique critical point is the origin (0, 0) which is a maximum 



(b) (3 > ^ for Xm<0,or^<f3<Y- for Xm > .■ three critical point (0, 0) 
and ±(x, —x) where x is a solution of: 



The origin is an inflection point while the other two are maximum points. 

(c) Y~ < f3 < -^j^ (Xm > 0): five critical points (0,0), ±{x,—x) and ±{x,x) 
where x is a solution of ^ and x is a solution of: 



point. 




(9) 




(10) 
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and $ is defined as a solution of: 



^"^ rill 



The origin is a minimum point, ±{x, —x) are maximum points and ±{x, x) are 
inflection points. 

(d) P > -^j^ (Xm > 0): nine critical points, ±{x,—x), ±{xl,X2), ±{x,x), ±{x2,xl), 
(0, 0) where x is a solution of x is a solution of / fi^) while for the other 
four points there is no analytical definition. The origin is a minimum point, 
±{x,—x) and ±{x,x) are maximum points. The other four can not be maxi- 
mum points. 

2. if Jii < and J12 > | Jii| (\m = Jii + J12 > 0; A™ = Jn - J12 < Oj 

(a) 0</3<i^.' the unique critical points is the origin (0, 0) which is a maximum 
point. 

(h) (3 > ^ ; three critical points (0,0) and ±(x,x) where x is a solution of 
The origin is an inflection point, ±(x,x) are maximum points. 

3. if Jn <0 andO < J12 < | Jii| (Xm = Jii + J12 < 0; A^ = Jn - J12 < Oj 
the unique critical points is the origin (0, 0) which is a maximum point. 

Corollary 1. Consider the Hamiltonian defined in (CP whit a = 1/2, Jn = J22 and 
hi = h2 = 0. If Jn and J12 have the same signs, the pressure functional f given by ^ 
admits a structure of critical points symmetric with respect to the vertical axis to those 
presented in Proposition 1. 

The corollary follows from Proposition 2 and the properties of symmetry previously 
exposed. 
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Proposition 3. Consider the Hamiltonian defined in (QP whit a = 1/2, Jn = J22 > 
and J12 < 0. 

1. As hi = h2 = the pressure of the system is given by: 



where f is given by x is a solution of ^ and t = 2(/3Aa/) ^ whit \m defined as 
in the previous proposition. 

2. The application of an infinitesimal field h = (/ii,/i2) with hi 7^ h2 selects the state 
such that the scalar product between the state and the field is positive; a field with 
equal components retains the symmetry of the system and does not allow to select a 
preferred status. 

For all other possible combinations of the parameters Ju and J12 similar results can be 
obtained simply by considering the model's symmetry, the Proposition 2 and its Corollary. 

We observe that as the external field is away and the bigger eigenvalue of the reduced 
interaction matrix is positive, at the temperature t = 1 the system undergoes a transition 
phase from zero to non-zero magnetization. 



Proof of Proposition 1 

We compute the exact solution of the thermodynamic limit exploiting a tail estimation 
on the number of configurations that share the same vector of the magnetization. In this 
way we obtain a lower and an upper bound for the partition function that converge to 
a same value as — )■ 00. This technique is used by Talagrand in [9j to compute the 
thermodynamic limit for the Curie- Weiss model. 
Denoted with <t; the configuration of the spins of the set Pi and: 




when t > 1 



when t < 1 



2 Proofs 




(12) 
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we can write: 

ZAr(J, h) = ^ exp {iSNgifii, /is)) . 

Lemma 1. Consider the set Q^i = 1}^' of all possible configuration cri. Let A^^ be 
a positive number defined by [W\). Then the following inequality holds: 



1 2^' 

--^exp(-iV,^(/.0) < A,^ < 2^'exp(-iV,^(/iO) (13) 



where C is a constant and is given by 

See appendix for the proof. This lemma allows to bound the partition function in the 
following way: 

--==^exp [[iNg{^) - N,,f{^{) - N2^{^^^)) < Z^(J,h) < 

< 2^ J^exp [pNgiti] - N,J^{fi,) - N2^{fi2)) 

then, since a sum of positive elements is always greater than its addends and always less 
than the largest addend times the number of elements, we have: 

1 2^ / \ / 

-^==exp iVmax/(M) < ^jv(J,h) < 2^(iVi + l)(iV2 + l)exp iVmax/(^; 

where the function / is defined in Hence for the pressure (jS]) we have: 
ln2 - ^( lnC + llnA^iA/'a ) +max/(/x) <pjv(J,h) < 



N\ 2 



< ln2 + 1 ln((iVi + 1){N2 + 1)) + max/(^). 

iV M 



Taking the limit for N ^ oo the Proposition 1 is proved. 
Proof of Proposition 2 

To study the space phase of the model defined by the Hamiltonian ([1]) under the hypothesis 
that a = 1/2, Ju = J22, hi = h2 = and the signs of Ju and Ju are opposite we 
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rescale the model's parameter with respect to the maximal eigenvalue Am of the reduced 
interaction matrix J in the following way: 

""'■^TW ^'-^in'' ^'-^iHin- 

\^m\ I^a/I pi Am I 

In this way the system of mean field equations ([8]) becomes: 




= tanh [j(aa;i + bx2) 
= tanh [j(&a;i + 0x2)] . 

By inverting the hyperbolic tangent in the two equations, we can write function 
of X2, and vice- versa X2 as a function of Xi. Therefore, when b ^ (that is the model 
doesn't degenerate toward two Curie- Weiss models) we can rewrite the equations in the 
following fashion: 

X2 = ■^(— axi -|- t tanh~^(a;i)) 
Xi = ^{—ax2 + ttaiah~^{x2))- 
Such a system lends itself to a graphic resolution. In fact the points that verify simultane- 
ously the two equations are the intersections of curves 71 and 72, respectively the graphs 
of the functions X2{xi) and Xi{x2)] 72 is the symmetric curve of 71 with respect to both 
the bisector of the Cartesian plane. In particular we will plot the curve 71 by studying the 
function X2{xi) and the curve 72 by symmetry. Before starting to analyze the different 
cases as the interaction parameters and the temperature change, we observe that X2{xi) 
is an odd function in the interval (—1, 1), it diverges as xi — )■ ±1 depending on the sign 
of b and its derivative are: 



= ( 1 1^2 - ' 4(0) = ^(t - a); (15) 



/// \ 2t Xi 

We are now ready to analyze in detail the different cases. 
1. Jii > and J12 < 

The maximal eigenvalue of the reduced interaction matrix is Am = Jii — J12, thus from 
(HM we have 0<a<l, — 1<6<0 and a — b = 1. In particular since b is negative. 



the function X2{xi) tends to +00 as xi —1 and to —00 as xi — )■ 1 and the sign of its 
second derivative is always opposite to that of the variable. 

(a)t>l (i.e. 0<f3<^) 

from (fT5|) we get that the function X2{xi) is monotonic decreasing and the value of its 
first derivatives in the origin is minor than —1. Thus the curve 71 lies above the bisector 
of the second and fourth quadrant when xi < and over when xi > 0, the symmetric 
curve 72 lies over the bisector when xi < and above when xi > 0, therefore the unique 
intersection is the origin. 



In order to understand what type of critical point is the origin, we consider the function 
/ defined in ([6]), we rescaled it following f lT^ : 



-1 




Figure 1: t=l.l, b=-0.8. The blue curve is 71, while the red is 72. 



/(Xi, ^2) = 2 ( ^ (I (^1 + ^2) + ^3:1X2^ - J^(Xi) - J^(X2) 



we compute its second derivatives: 
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and the determinant of its Hessian matrix: 



a-^—^]~b']. (16) 



In the origin we have: 



Hf{0, 0) = ^((« - -b') = ^,iib + l- tf - 6^) = 
d'^f 1 

Thus, since t>l,0<a<l and — 1 < 6 < the origin is a maximum point of /. 
(b)t = l (t.e. (5 = ^) 

the function X2{xi) is still monotonic decreasing and the value of its first derivatives in 
the origin is equal to —1. Therefore the position of the curves 71 and 72 with respect to 
the bisector of the second and fourth quadrant are the same of previous case. And so 
there is only one intersection between them in the origin. The determinant of the Hessian 
matrix ( fT6l) in this point is equal to zero. To understand what type of critical point 
is the origin we consider the following change of variable that diagonalizes the reduced 
interaction matrix: 

_ X1+X2 

(17) 




2 

The Taylor expansion of the function / expressed in the new variable assures that such a 
point is a maximum point. (For the details see the Appendix). 

{c)a<t<l (^.e. < P < ^-J 

the function X2{xi) is still decreasing but its derivative in the origin f lTSj) is bigger than 
—1, thus the curve 71 intersects the bisector of the second and fourth quadrant (and 
therefore the curve 72) in (0, 0) and in other two points symmetric with respect to the 
origin, (— a;,x) and (x, — x), where x is the positive solution of the equation v{x) = 0: 

f (x) = ttanh~"^(x) — X. (18) 
11 



Figure 2: t=0.75, b=-0.8. The blue curve is 71, while the red is 72. 



Since the function, v given by f llSp . intersects the horizontal axis in the origin, reaches its 
minimum in x = a/1 — t and goes to infinity for large x, it follows: 



X > VT^. (19) 

This inequality turns to be useful to evaluate the nature of the critical points ±(x, —x). 
In these points the determinant of the Hessian matrix of the function / f|T6|) is: 

Hf{-i,S:) = Hf{i,-i) = ^(^(^a~^-^^ - b"^. (20) 

Since in the considered range of temperature the function on the r.h.s of fl20l) is monotonic 
increasing for positive x, by the inequality f|T9|) we have: 

Hf{-x,x) = Hf{x,-x) > ^{{^- 

On the other hand the first term of the Hessian matrix in ±( decreasing function 

and thus: 

-X, X = x, —X = — a — -r < 



dx"^ ' dx"^ ' 2t \ 1 — x^ 

< — a - 1 < 0. 

2V ' 
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Then ±(x, —x) are maximum points. On the contrary the origin is an inflection point 
because the determinant of the Hessian matrix ( IT6|) computed in such a point is negative 
for the considered range of temperature. 

(d) 2a - 1 < t < a with 2a - 1 > or < t < a (i.e. > and -j^ < P < ^ or 
A„ < and (3 > j^) 

The function X2{xi) has a minimum point in Xm = — a/1 — t/a and a maximum point 
in Xm = ^/l — t/a. Therefore the curves 71 and 72 intersect each other in at least three 
points: the origin and ±(x, —x), where x is again the positive solution of v{x) = 0, with 
V defined in flTSl) . However, since in the origin the first derivative of X2{xi) is positive and 
less than one, 71 lies above the bisector of the first and third quadrant when xi < and 
over when xi > 0. Instead 72 lies over the bisector when xi < and above when xi > 0. 
This scenario excludes the possibility of additional intersections between the two curves. 





1- 






\ 0.6- 
\ 0.4- 
\ 0.2- 




i -oS -0^6 -0^4 -a.z'yV 
/ """"^ 

/ -o>-. 


0!2 QA 0.6 0^^ 


1 


/ -1- 







Figure 3: t=0.5, b=-0.3. The blue curve is 71, while the red is 72. 

Again ±(x, —x) are maximum points and the origin is an inflection point, 
(e) t = 2a - 1 with 2a - 1 > (i.e. > and = j^) 

in this case the only difference from the previous one is the value of the flrst derivative of 
the function X2{xi) at the origin, which now is exactly equal to 1. This does not imply any 
qualitative change in the relative positions of curves 71 and 72, that continue to intersect in 
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the points {—x, x), (0, 0) and (x, —x). The analysis of the Hessian matrix of the function 
/ assures that ±(x, —x) are still maximum points but does not give information about 
the nature of the origin. Considering again the change of variable fll7l) we can conclude 
that the origin is an inflection point. (See the Appendix). 

(f) t < t < 2a — 1 tvhere t is a solution of ttanh^^{Vl — i) — — i = 1h\/\ — t 

(i.e. Am > and < /3 < -57— where B is a solution of 071) 
the function X2{xi) has the same intervals of monotony of the previous cases but the 
value of its first derivatives in the origin is bigger than 1. This implies that it intersect 
the bisector of the first and third quadrant in two other points beyond the origin, (— x, —x) 
and where x is the positive solution of the equation v{x) = 2bx, where v is the 

function defined in (fT8|) . Thus the curves 71 and 72 intersect each other in at least five 
points: (0, 0), ±(x, —x), ±(x, x). The graphic representation presented in [2] excludes the 
possibility of further solutions that will appear as t < t. 



The computations made on the Hessian matrix assures that the origin is a minimum 
point and ±(5;, —x) are maximum points. The new solutions ±{x,x) are infiection points 
because in these points the determinant of the Hessian matrix of the function / can be 



-1 




Figure 4: t=0.3, b=-0.3. The blue curve is 71, while the red is 72. 
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expressed in the following way: 

2 




6' 



2 



= ^ - l) 

where — 1 < {x2{x)) < 1 in the considered range of temperature. To prove this statement 
we consider the inequality: 

x> \ll ^ 



1 + 26 

true for the solutions of the equation v{x) = 2bx for the same reason for which the 
inequality ( fT9l) is true for the solution of v{x) = 0. Since: 



2b 



and the function x'2{xi) is monotonic decreasing for positive values of the abscissa, we have 
that X2{x) < 1. In particular the value of this derivative, initially positive, decreases as the 
temperature t decreases till the values —1 that it reaches for t = i. At this temperature 
the determinant of the Hessian matrix computed in ±(x, x) is equal to zero. Nevertheless, 
without need to resort to the diagonalization of the reduced interaction matrix, we can 
say that these points, equal to ±(-\/l — t, y/1 — t), are inflection points. In fact: 



^(±v^, ±VT^t) = ^(±v^, ±VT^t) = — ^ ^ 

and the other partial derivatives of the third order of the function / are equal to zero, 
(g) < t < t (i.e. A„ > and (3 > ) 

the graphical representation presented in p] shows that as the temperature decreases 
below the value i a further solution {x^, X2) of the mean field equations appear, with 
X2>x\> 0; the other three solutions are generated by the symmetry with respect to 
both the bisector. In these points it is not always possible to detect the sign of the 
determinant of the Hessian matrix. In any case, since they are linked to a maximum 
point by a monotonic piece of one of the two curves 71 and 72, we can exclude that they 
are in turn maximum points. (For further details see [5]). When t is small enough it can 
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be shown analytically that they are inflection points. 
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Figure 5: t=0.25, b=-0.3. The blue curve is 71, while the red is 72. 

As in the previous case the origin is a minimum point and ±(a:, —x) are maximum points; 
on the contrary ±(a;,a;) became maximum points. In fact since x'^ix) < — 1 as t < t, the 
determinant of the Hessian matrix of / in these points is positive and: 

2. Jii < and J12 > |Jii| 

The maximal eigenvalue of the reduced interaction matrix is Aa^ = Ju + Ju, thus from f|T^ 
we have a < 0, 6 > 1 and a + b = 1. In particular since b is positive, the function X2{xi) 
tends to —00 as Xi — )■ —1 and to -|-oo as a;i — > 1 and the sign of the second derivative of 
X2{xi) is always the same of the variable. Moreover the signs of a and b assure that the 
function X2{xi) is always monotonic increasing. Therefore we can distinguish only two 
scenarios: 

(a)t>l r^.e. 0<P<^) 

the value of the first derivative of X2{xi) computed in the origin is major or equal to 1. 
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Thus 7i lies over the bisector of the first and third quadrant when xi < and above when 
Xi > 0. Therefore since 72 is symmetric to 71 with respect to this bisector the curves 
intersect each other once in the origin. Analyzing the Hessian matrix of the function / in 
this point for t > 1 we have: 



Thus (0, 0) is a maximum point. 

As t = 1 in the origin, the determinant of the Hessian matrix of the function / is equal 
to zero, however considering the change of variable f[T7|) and the Taylor expansion of the 
function / we get that such a point is a maximum point. (See the Appendix). 

(b)0<t<l (z.e. (3>^) 

the value of the first derivative of X2{xi) in the origin is less than 1 and then the curves 
7i and 72 intersect each other in three point: (0, 0) and other two points symmetric with 
respect to the origin (— x, — x), {x,x) where x is the positive solution of the equation 
v{x) = 0. 

By inequalities f lT9|) we have: 




1 

1 , 



(a-t) < 0. 



(t-l)(t-l + 26) > 



2t 




Hf{—x, —x) = Hf{x, x) 




) 



Therefore ±(x, —x) are maximum points. On the other hand since: 



1 

4^2 



(t-l)(t-l + 26) <0 



the origin is an inflection point. 
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3. Jii < and < J12 < |Jii| 

The maximal eigenvalue of the reduced interaction matrix, Am = Jn + J12, is negative. 
Therefore by (IT^ we have a < — 1, 6 > and a + b = —1. Since the signs of parameters 
a and b are the same of the previous case, the comments made about the limit behavior, 
the monotony and the concavity of the function X2{xi) are still valid. Furthermore, since 
a > —1 the value of the first derivative of X2{xi) computed in the origin is bigger than 1 
for all if: > 0. As we have seen before this imply that curves 71 and 72 intersect each other 
only in the origin. Since: 

Hfi0,0) = ^it+l)it + 2b+l) > 
g^(0,0) = -(a-t) <0, 
such a point is a maximum point. 



Proof of Proposition 3 

1. If 2a — 1 < the statement follows directly from Proposition 2. If 2a— 1 > the same 
is true only for t > i. For < t < t the function / admits four maximum points: 
±(5;, —x) and ±(x, x), therefore to find what is the absolute maximum between them 
we have to compare the values of the function in these points. For this purpose it 
is more convenient to consider the function P{xi,X2) = tf{xi,X2)'. 

P{xi, X2) = ^a{xl + xl) + ^6x1X2 - ^ (J^(a;i) + ^{x2)) ■ 

Since the function P is pair we choose to compute the value it assumes in the points 
with positive abscissa, (x, —x) and (x, x) where x and x are positive: 

P(x, — x) = -ax^ — 2^^^ ~ tJ^(x) = -x^ — t^{x) 
P(x, x) = ^ax^ + ^6x^ - tJ^(x) = ^(a + b)x^ - tJ{x). 

We observe that in the points (x, — x) and (x, x), the first and the second derivatives 
with respect to the temperature of the function P can be write: 

dP(M(t)) ^, , d^PiiMit)) ^, , , , ^ 
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where M{t) denotes in a generic way one of the two maximum points and XM{t) 
is its abscissa. Since the function is pair, convex and monotonic increasing for 
positive abscissa, and xm(^) is decreasing, the function P is decreasing ad concave 
along the maximun points. As the temperature tends to zero both x and x tend to 
one, therefore: 

, ,1 1 - -\ 2a — 1 

hmPix, —x) — limPia;, a;) = — - — . 

Since < a < 1 we obtain that as t — > the pressure reaches the higher value in 
(x, —x). Also as i = i it is easy to show that: 

P{x,x) < P{x, -x). 

Thus, since x{t) > x{t) as < t < t we conclude that: 

P{x, -X, t) > P{x, X, t) for < i < t. 

2. If we apply an infinitesimal field h = (/ii, /i2) to the considered system, the maximum 
points of the function / remain around the points that were maximum points when 
the field was away. We are interested in showing the field's ability to select a 
preferred state, thus we omit the case t > 1 in which there is only one possible 
state. On the contrary for other temperature values the function / admits at least 
two maximum points. The gap between the values that the function P assumes in 
the points ±(x, —x) and ±{x,x) in absence of the field assures that when the field 
is infinitesimal the maximum points close to ±{x,x) cannot be the effective state 
of the system. In order to understand which of the two remaining points is the 
preferred we distinguish three cases. 

If the field components are equal, h = {hi, hi), the mean field equations and their 
solutions retain the symmetry with respect to the bisector of the first and third 
quadrant that they presented as the field was away. Therefore for each temperature's 
value the pressure computed in the two global maximum points is the same and the 
field h is not able to select a preferred state between them. 
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If the field components arc opposite, h = {hi, —hi), the system retain the symmetry 
with respect to the bisector of the second and fourth quadrant, thus we have two 
solutions s^^^ and s^^^ on this bisector; their first coordinates are respectively the 
positive and negative solution of the mean field equation: 

Xi — tanh ^- {xi + hi) 

Therefore between s^^^ and s^'^\ that are on opposite side with respect to the origin, 
the field selects the point whose direction is identical to its own. 
At least we consider the case in which the field components are generic, that is 
h = (/ii,/i2) with \hi\ 7^ |/i2|. Decomposed the field along the eigenvector of the 
reduced interaction matrix J: 

{hi,h2) = {h',-h') + {h",h") 

we apply the field components in two steps. Firstly considering only {h', —h'), as 
we have seen above, we obtain two solutions s^^^ = {x^^\ —x^^^) and s^"^^ respectively 
close to {x,—x) and {—x,x). The subsequent application of {h",h") causes an 
infinitesimal shift of s^^^ and s^^^ in this direction. We look for the coordinates of 
s^^^ in the form: 

x?(t)^x^'\t) + x';(t) 
41) (t) = -x(^){t) + x';{t) 

where x'l is the infinitesimal shift in the diagonal direction. Prom the mean field 
equation we obtain: 

xW + x'/ = tanh [| (x« + (a + + /i")] 
-xW + x'l = tanh [i (-x(i) + (a + b)x'l + h")] . 

Developing the hyperbolic tangent and neglecting the infinitesimal terms of second 
order of the type x"h" we obtain two equations both equal to the following: 

x'l = (1 - tanh 



a + b 



-x'{ + 



h" 
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whose principal solution is those that has the same sign of h" . Since the same 
argument applies to the second solution we conclude that the field h = {hi, with 
I /ill 7^ |/i2| moves ±(x, — x) in two solutions which remain on opposite half-plane 
with respect to the bisector of the first and third quadrant. Furthermore the field 
ability to broke the symmetry of the function / ensures that the absolute maximum 
point is the one whose scalar product with the field is positive. 

3 Comments 

In this paper we have deeply investigated the behaviour of the symmetric bipartite 
mean-field model as the external field is away or small. In particular, by studying the 
nature of the critical points of the pressure functional we have highlighted for which values 
of parameters the model undergoes a phase transition. In addition to being interesting 
in themselves, the obtained results can provide good tools to describe more complex 
models, like the bipartite non-symmetric mean-field model, already used in the statistical 
approach to socio-economical sciences [3], or the multipartite symmetric mean- field model, 
that seems to have a very interesting behaviour with non trivial phase transitions. These 
models will be subject of further investigation of the authors. By analyzing them is 
expected to relate the critical and non-critical phases to different limiting behaviours of 
the sums of the spins of such models [1]. 

Appendix 

A Proof of lemma [1] 

As mi{(T) = fii, the configuration cr/ contains A^/(l + fii)/2 times the value 1 and 
Ni{l — fii)/2 times the value —1, thus we have: 
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Using Stirling's formula, n\ ~ n"e "\/27m, we get: 

^ [2 1 ^ 



1 2^' 1 

> 



1 2^' 



In this way we obtain a lower bound of ^4^^ . 

To obtain an upper bound for A|^^ we suppose the spins Ui are independent Bernoulli 
random variable. In this case P{ai = 1) = P{cFi — —1) = 1/2 and all configurations ai 
have same probability, hence: 

where by definition of the magnetization: 



Take A > by Chebyshev's inequality we can bound the above probabihty: 

> f^iNi \ < e-^«^' n£;[exp(Aa,)] = 

— exp{Ni{—XijLi + In cosh A)) < 

< mm{exp{Ni{— Xni + In cosh A))} 

A 

where E denotes the expectation value. 

If I//.; I < 1, the previous exponent is minimized for: 

-1/ X 1. fl+fil 



A = tanh-^(/.,) = -ln(^^). (21) 



Since 1 / (cosh^ y) = 1 — tanh^ y the following equahty holds 

1 



lncoshA = --ln(l-/xf). (22) 



22 



Thus by (|2T1) and (122]): 

mm{exp{Ni{—Xfii + In cosh A))} = exp{—NiJ^{fj,i)). 

A 

Hence the upper bound for is: 

B Diagonalization 

In the proof of proposition 1 we have seen that in three cases the determinant of the 
Hessian matrix of the function /, given by (|6]), is equal to zero when computed in the 
critical point (0,0). In particular this happen as Jn > and Ju < for t = 1 and 
t = 2b + 1 and as Jn < and J12 > |Jii| for t = 1. Thus, in order to determinate 
the nature of the origin we consider the change of variable that diagonalize the reduced 
interaction matrix J: 



and we compute the partial derivatives of / with respect to the new variables. Obviously 
the derivatives of the first order: 




X-Y 



X + Y 



we apply it to the function / 




dX 
dl 
dY 



iX,Y) 



{X,Y) 



a + b 



t 

a — b 



t 



X - ^ (tanh-^(X -Y)+ tanh-^(X + Y)) 
■Y -- (tanh"^(X + F) - tanh"^(X - Y)) 



are equal to zero in the origin. The derivatives of the second order are: 
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In particular in (0, 0) we have: 

d^f ,^ a + b 



' ^ t 

■' (0,0) = 



dXdY 

Thus, since a — 6 = 1 as Jn > and Ju < and a + 6 = lasJii<0 and Ju > \Jii\, 
in all the three cases the determinant of the Hessian matrix is equal to zero. This is the 
reason why we need to compute partial derivatives of order higher than the second. Those 
of the third order are: 



dX^^ ' ' dXdV^^ ' ' 2\{1-{X + Y)Y (l-(X-F) 
9'f(XY)-^^(XY)- V 2(^ + ^) I 



2\2 



dY^' ' ' dX^dY^ ' ' 2 {X + Y)^)^ {1-{X-Y)^)\ 

It is easy to check that these derivatives computed in (0,0) are equal to zero. At least 
the partial derivative of the fourth order are: 

Q4j Q^j g4j / I 

1 4(x + y)2 4(x-y)2 



;i-(X-y)2)2 (1 _ (X + y)2)3 (i-(x-y)2)3 



a4 ^ q4 X 



4(x + y)2 4(x-y)^ 



+ 



(i-(x-y)2)2 (i_(x + y)2)3 (i-(x-y)2)3 

In particular in (0,0) these derivatives become: 



(0,0) = ^(0,0) = ^^^(0,0) = -2 

0. 



V ' ^ 9y4^ ' ^ 9X29y2 



dX^dY dXdY^ 

Now we can approximate the function / in a neighborhood of the origin with its Taylor 
expansion till the fourth order. In particular as Jn > and J12 < and i = 1 we have: 

f{X, Y) = 2bX^ -^{X + Yf + o{X^ + y"). 
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Since for these values of the parameters b is negative, this polynomial is negative. There- 
fore the origin is a maximum point. As Jn > and J12 < and t — 2b + 1 the Taylor 
approximation of / is: 

f{X, Y) = -2bY^ -L(^x + Yf + o{X^ + Y^). 

It is easy to check that this polynomial is negative when computed in the point of the 
horizontal axis and positive in the point of the vertical axis. Thus the origin is an inflection 
point. At least as Jn < and J12 > | Jii| and t — 1 we have: 

f{x, Y) = -2bx^ - ^(x + y)' + o(x' + y"). 

In this case b is positive, thus this polynomial is negative and the origin is a maximum 
point. 
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